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We introduce a semiclassical Einstein-Langevin equation as a consistent dynamical equation for a 
first order perturbative correction to semiclassical gravity. This equation includes the lowest order 
quantum stress-energy fluctuations of matter fields as a source of classical stochastic fluctuations 
of the gravitational field. The Einstein-Langevin equation is explicitly solved around one of the 
simplest solutions of semiclassical gravity: Minkowski spacetime with a conformal scalar field in its 
vacuum state. We compute the two-point correlation function of the linearized Einstein tensor. This 
calculation illustrates the possibility of obtaining some "non-perturbative" behavior for the induced 
gravitational fluctuations that cannot be obtained in perturbative quantum gravity. 
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Semiclassical gravity describes the interaction of the 
gravitational field, which is treated classically, with quan- 
tum matter fields. The theory is mathematically well de- 
fined and fairly well understood, at least for linear matter 
fields . The equation of motion for the classical met- 
ric is the semiclassical Einstein equation, which gives the 
back reaction of the matter fields on the spacetime; it is a 
generalization of the Einstein equation where the source 
is the expectation value in some quantum state of the 
matter stress-energy tensor operator. 

In the absence of a complete quantum theory of grav- 
ity interacting with matter fields from which the semi- 
classical theory can be derived, the scope and limits of 
semiclassical gravity are less well understood. It seems 
clear, however, that it should not be valid unless gravi- 
tational fluctuations arc ncglig ibly small §{§. This con- 
dition may break down when the matter stress-energy 
has appreciable quantum fluctuations ppj, given that 
a quantum metric operator should couple to the stress- 
energy operator of matter and, thus, fluctuations in the 
stress-energy of matter would induce gravitational fluc- 
tuations In recent years, a number of examples have 
been studied, including some quantum fields in cosmolog- 
ical models and even in flat spacetimes with non-trivial 
topology, where, for some states of the fields, the stress- 
energy tensor have significant fluctuations [||. It thus 
seems of interest to try to generalize the semiclassical 
theory to account for such fluctuations. 

In this paper we propose a generalization of semiclas- 
sical gravity based on a semiclassical Einstein-Langevin 
equation, which describes the back reaction on the space- 
time metric of the lowest order stress-energy fluctua- 
tions. This results in an effective theory which predicts 
linear stochastic corrections to the semiclassical metric 
and may be applicable when gravitational fluctuations of 
genuine quantum nature can be neglected. It should be 
stressed that the gravitational fluctuations predicted by 
this stochastic semiclassical theory are "passive" rather 
than "active", i.e., they are induced by the matter field 



stress-energy fluctuations [Q. Thus, we cannot expect 
that stochastic semiclassical gravity gives a satisfactory 
description of gravitational fluctuations in all situations. 
In spite of that, this theory may have a number of in- 
teresting applications in the physics of the early universe 
and of black holes, where one may expect significant mat- 
ter stress-energy fluctuations. 

The idea which motivates such stochastic theory is 
that, in the transition from the full quantum regime to 
the semiclassical one, there should be some mechanism 
which effectively suppresses the quantum interference ef- 
fects in the gravitational field. Once such decoherence 
process has taken place, gravity might undergo an inter- 
mediate regime in which it would fluctuate, but it would 
do so classically. If such a regime exists, an effective 
probabilistic description of the gravitational field by a 
stochastic metric field may be possible Q. 

Another motivation for our work is to connect with 
some recent results on equations of the Langevin-type 
that have appeared in the context of semiclassical cosmol- 
ogy and which predict stochasti c per turbations around 
some cosmological backgrounds [§,0. Although phys- 
ically motivated |l]J| , the derivation of these equations, 
obtained by functional methods, is formal and doubts 
may be raised on the physical significance of the pre- 
dicted fluctuations. Our approach is complementary to 
these functional methods, it allows to write these equa- 
tions in a general form and it links the source of such 
metric fluctuations with the matter stress-energy fluctu- 
ations. 

The idea, when implemented in the framework of a per- 
turbative approach around semiclassical gravity, is quite 
simple. We start realizing that, for a given solution of 
semiclassical gravity, one can associate the lowest order 
matter stress-energy fluctuations to a classical stochastic 
tensor field. Then, wc seek a consistent equation which 
incorporates this stochastic tensor as a source of a first 
order correction to semiclassical gravity. It is important 
to remark that, even if this approach differs from those 
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based on functional methods, the semiclassical Einstein- 
Langevin equation introduced here can actually be for- 
mally derived using those methods. The details of such 
derivation will be given in a subsequent paper 

As a simple application of stochastic semiclassical grav- 
ity, we have explicitly solved the Einstein-Langevin equa- 
tion around a solution of semiclassical gravity consisting 
on Minkowski spacetime with a conformal scalar field in 
the Minkowskian vacuum state. We have computed the 
two-point correlation function for the induced linearized 
Einstein tensor. Even if, as expected, such correlation 
function has negligible values for points separated by 
scales larger than the Planck scales, it is interesting to 
compare our result with the one obtained from a pertur- 
bative approach to quantum gravity |ll| . The result hints 
at the possibility of obtaining some non-perturbative (in 
the Planck length) behavior for the induced gravitational 
fluctuations which could not be obtained in perturbative 
quantum gravity. Throughout we use the (+ + +) sign 
conventions and work in units in which c = fi = 1. 

The semiclassical Einstein-Langevin equation. — We 
start with the semiclassical Einstein equation. Let 
(M, g ab ) be a globally hyperbolic four-dimensional space- 
time and consider a linear quantum field $ on this back- 
ground. Working in the Heisenberg picture, let and 
p[g] be respectively the field operator and the density op- 
erator describing the state of the field. Such state will 
be assumed to be physically acceptable in the sense of 
Ref. 0. 

The set (M, g ab , <1>, p) is a solution of semiclassical grav- 
ity if it satisfies the semiclassical Einstein equation: 

^ (G ab [g] + Ag ab ) - 2 (aA ab + f3B ab ) [g] = (f£)[g], 

(1) 

where T^[g\ is the renormalized stress-energy tensor op- 
erator for the field $ [g] , which satisfies the corresponding 
field equation on the spacetime (M,g ab ) 7 and the expec- 
tation value is taken in the state described by p[g] 
In the above equation, 1/G, A/G, a and (3 are renormal- 
ized coupling constants, G ab is the Einstein tensor, and 
A ab and B ab are the local curvature tensors obtained by 
functional derivation with respect to the metric of the 
action terms corresponding to the Lagrangian densities 
C abc dC abcd and R 2 , respectively, where C abc d is the Weyl 
tensor and R is the scalar curvature. A classical stress- 
energy tensor can also be added to the right hand side of 
Eq. (|l|), but, for simplicity, we shall ignore this term. 

Given a solution of semiclassical gravity, the matter 
stress-energy tensor will in general have quantum fluc- 
tuations. To lowest order, such fluctuations may be de- 
scribed by the following bi-tensor, which we call noise 
kernel. 

8N abcd (x,y) = ({i ab {x), t cd (y)})[g], (2) 



where { , } means the anticommutator and i ab = T ab — 
(T ab ), where denotes the unrenormalized stress- 

energy "operator" [we use the word "operator" for T ab in 
a formal sense; it should be understood that the matrix 
elements of this "operator" are suitably regularized and 
that the regularization is removed after computing the 
right hand side of (Q)]. For a linear matter field, this 
noise kernel is free of ultraviolet divergencies, and the 
"operator" i ab in (||) can be replaced by the operator 

We want now to introduce an equation in which the 
stress-energy fluctuations described by (||) are the source 
of classical gravitational fluctuations, as a perturbative 
correction to semiclassical gravity. Thus, we assume that 
the gravitational field is described by g ab + h ab , where 
h ab is a linear perturbation to the metric g ab , solution 
of Eq. (Q) . The renormalized stress-energy operator and 
the density operator describing the state of the field will 
be denoted as T^g + h] and p[g + h], respectively, and 

T^)[g + h] will be the corresponding expectation value. 
In order to write an equation which describes the dy- 
namics of the metric perturbation h abl let us introduce a 
Gaussian stochastic tensor field t; ab characterized by the 
following correlators: 

(U(x)) c = 0, (U{x)U(y)) c = N abcd (x,y), (3) 

where ( ) c means statistical average. Note that the two- 
point correlation function of a stochastic tensor field £ a b 
must be a symmetric positive semi-definite real bi-tensor 
field (since, obviously, (£ ab (x)£ cd (y)) c = (£cd(y)£ab(x)) c )- 
Since T5 is self-adjoint, it is easy to see from the defi- 
nition (|2j) that N abcd (x, y) satisfies all these conditions. 
Therefore, relations (g), with the cumulants of higher or- 
der being zero, do truly characterize a stochastic tensor 
field £ a fc. One could also seek higher order corrections 
which would take into account higher order stress-energy 
fluctuations, but we stick, for simplicity, to the lowest 
order. The simplest equation which can incorporate in 
a consistent way the stress-energy fluctuations described 
by N abcd (x, y) as the source of metric fluctuations is 

(G ab [g + h] + A (g ab + h ab )j 

- 2 (aA ab + /3B ab ) [g + h] = (f£) [g + h] + 2£ ab , (4) 

which must be understood to linear order in h ab . This is 
the semiclassical Einstein-Langevin equation, which gives 
a first order correction to semiclassical gravity. Notice 
that, in writing Eq. (Q), we are implicitly assuming that 
h ab is also a stochastic tensor field. 

We must now ensure that this equation is consistent 
and, thus, it can truly describe the dynamics of metric 
perturbations. Note that the term £ ab in Eq. (||) is not 
dynamical, i.e., it does not depend on h abl since it is de- 
fined through the semiclassical metric g ab by the correla- 
tors (0). Being the source of the metric perturbation h ab , 
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this term is of first order in perturbation theory around 
semiclassical gravity. Let us now see that £ a j, is covari- 
antly conserved up to first order in perturbation theory, 
in the sense that the stochastic vector field V a £ab is de- 
terministic and represents with certainty the zero vector 
field on M (y a means the covariant derivative associ- 
ated to the metric g a b)- Using that X7 a T a F {g] — 0, we 
have \/ a N a bcd{x,y) = and, from the covariant deriva- 
tive of the correlators (||), we get (v a £a&)c = and 
(V^ab(x)\7^cd(y))c = 0. It is thus consistent to in- 
clude the term in the right hand side of Eq. (Q). 
Note that for a conformal matter field, i.e., a field whose 
classical action is conformally invariant, the stochastic 
source £ a b is "traceless" up to first order in perturba- 
tion theory. That is, g ab £, a b is deterministic and repre- 
sents with certainty the zero scalar field on M. In fact, 
from the trace anomaly result, which states that g ab T^\g] 
is in this case a local c-number functional of g c d times 
the identity operator, we have that g ab (x)N a bcd(x,y) — 
0. It then follows from (g) that (g ab ^ a b)c = and 
(g ab (x)£,ab(x)g cd (y)£,cd(y))c = 0. Hence, in the case of 
a conformal matter field, the stochastic source gives no 
correction to the trace anomaly. 

Since Eq. (^) gives a linear stochastic equation for 
h a b with an inhomogeneous term ^ a0l a solution can be 
formally written as a functional h ao [£] of the stochastic 
source £ cc j. Such a solution can be characterized by the 
whole family of its correlation functions. By taking the 
average of Eq. (pj), one sees that the metric g a b + (h a b) c , 
must be a solution of the semiclassical Einstein equation 
linearized around g ao . For the solutions of Eq. (Q) we 
have the gauge freedom h ab -> h' ab = h ab + VaC& + V&Ca, 
where C a is any stochastic vector field on M which is 
a functional of the Gaussian stochastic field £ C£ ;, and 
Ca = 3afcC & i so that h a b and h' ab are physically equiva- 
lent solutions. 

Vacuum fluctuations in fiat spacetime. — As an ex- 
ample, we shall now consider the class of trivial solutions 
of semiclassical gravity. Each of such solutions consists 
of Minkowski spacetime, (M 4 , rj a b), a linear matter field, 
and the usual Minkowskian vacuum state for this field, 
p[rj] = |0)(0|. As it is well known, we can always choose a 
renormalization scheme in which (0| T^ b |0) rj\ = 0. Thus, 
each of the above sets is a solution to Eq. (l|) with A = 0. 
Note that, although the vacuum |0) is an eigenstate of 
the total four-momentum operator which can be defined 
in Minkowski spacetime, such state is not an eigenstate 
oiT^[rj\. Hence, even in these trivial solutions, quantum 
fluctuations of the matter stress-energy are present, and 
the noise kernel defined in (g) does not vanish. This fact 
leads to consider the Einstein-Langevin correction to the 
trivial solutions of semiclassical gravity Jl2] . For a mass- 
less field, these vacuum stress-energy fluctuations should 
be characterized by the only scale available, the Planck 
length, and, thus, they should be very small on macro- 
scopic scales, unlike the cases considered in Refs. (6). 



The corresponding Einstein-Langevin equation be- 
comes simpler when the matter field is a massless con- 
formally coupled real scalar field. In the global inertial 
coordinate system {ie m }, the components of the flat met- 
ric are simply rj^v — diag(— 1, 1, 1, 1). We shall use Gfy, 
A^l and Bfy to denote, respectively, the components of 
the tensors G a t>, A a b and B ao linearized around the flat 
metric. These tensor components can be written in terms 
of C?'f' as 



7" 



aw — _(v nmot _ f:ar^m\ 



Bpv — avG 



(1)Q 

a J 



where J~pv is the differential operator T ^ = rj^ v 

v • 

It is also convenient to introduce the Fourier transform, 
f(p), of a field f(x) as f(x) ee (2^)- 4 fd 4 p ^ /»■ 
Eq. (0) reduces in this case to 
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G™(x)-2(aA$ + PB$)(x) 



+ J d 4 yH(x-y;^)A^(y) = 2^(x), 



(5) 



where a = a + 1/(3600tt 2 ), j3 = - 1/(34560tt 2 ), and 
H{p;H 2 ) ee (ln|p 2 /p 2 | -iTrsignp 6>H? 2 )) /(19207T 2 ), 
being fi a renormalization mass scale [ p.0| . The compo- 
nents £ M „ of the stochastic source tensor satisfy in this 
case 



(^{X)ial3(y))c 



~ FpvaP N(x - 



y), 



where Ty, va p = iT^T^ - T^ u T a p and N(p) = 
6»(-p 2 )/(19207r). 

Equations can be solved for the components G™ 
of the linearized Einstein tensor. Using a procedure sim- 
ilar to that described in the Appendix of Ref. || , we find 
the family of solutions which can be written as a lin- 
ear functional of the stochastic source and whose Fourier 
transform G l *l(p) depends locally on £ a p(p)' Each of 
such solutions is a Gaussian stochastic field and, thus, 
it can be completely characterized by its one-point and 
two-point correlation functions. The one-point correla- 
tion functions, i.e., the averages (Gj^) c , are solution 
of the linearized semiclassical Einstein equations ob- 
tained by averaging Eqs. (||); solutions to these equations 
were first found by Horowitz [jL4|. We can then com- 
pute the two-point correlation functions Guvap{x,x') = 
(G« 1 >(^)Gi 1 »(x')) c -(G< 1 )(x)) c (G^( a :'))c & These cor- 
relation functions are invariant under gauge transforma- 
tions of the metric perturbations and give a measure of 
the induced gravitational fluctuations in the present con- 
text, to the extent that these fluctuations can be de- 
scribed by exact solutions to Eqs. (||) (the result would 
be different if one applied some "reduction of order" pro- 
cedure H). We get 



Gnvap{x, x') 
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G 2 Fuvu0 Q{x — x), 



(6) 
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where Q(p) — 6(—p 2 ) 1 + 16nGp 2 H(p; [i 2 ) , with /x = 
[i exp(19207r 2 a). Introducing the function <^(x; A) 
I - xln(Ax/e)] 2 + ttV, with x > and A > 0, Q(x) 
can be written as 



Q(x) = 



(120tt) 3 / 2 1 



2tt 3 L3 



OO 

: / dq° cos 
/o 



d|q||q| sin 



Vl20^ 

r I ? 



^Hz 2 ;A)' 



(7) 



where x" = (x°,x) and ^ = (q°,q), A = 1207re/(Lp£ 2 ), 
and Lp = \/G is the Planck length. If we assume 
that p, < Lp 1 , then A > 10 3 . Let us consider the 
case when (x — x')^ is spacelike, then we can choose 
(x — x')^ = (0, x — x') and Q liVa p{x,x') will be a func- 
tion of x — x' only. From @ and p), one can see that 
5oooi(x — x') = £oijfc(x — x') = 0. The remaining compo- 
nents are non-null and can be explicitly computed per- 
forming the integrals with some approximations. The 
result depends on the constant k = ln(Axo(A)/e), be- 
ing Xo(A) the value of x where ip(x',X) has its mini- 
mum (this can be found by solving the equation 7r 2 xo 
I — xoln(Axo/e)] [1 + ln(Axo/e)] numerically). For x ^ 
x', we get 



j'OOOO 



(x- 



1 



4 
a 



12 

^2 



24 

^3 



24 

7^4 



where a = b |x — x'|/Lp, 
b = (4a/7T 2 )[l57r(V7^ 
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I- (2 /it) arctan(K/-7r) and 

7r 2 — k)] . We should remark 
that this is not an expansion in a. Similar results are 
obtained for the other non-null components. From these 
results, one can conclude that, as expected in this case, 
there are negligibly small correlations for the Einstein 
tensor at points separated by distances large compared to 
the Planck length. Thus, at such scales, the semiclassical 
approach is satisfactory enough to describe the dynam- 
ics of gravitational perturbations in Minkowski space- 
time 0J|. Deviations from semiclassical gravity start 
to be important at Planckian scales. At such scales, 
however, gravitational fluctuations of genuine quantum 
nature cannot be neglected and, thus, the classical de- 
scription based on the Einstein-Langevin equation would 
break down. Note, however, the factor e~ CT in our result, 
which is non-analytic in the Planck length and gives a 
characteristic correlation length of the order of Lp. This 
kind of behavior cannot be obtained from a perturba- 
tive approach to quantum gravity, in which one expands 
physical quantities as a power series in l? p |l3|]. Actu- 
ally, if we had naively assumed that (|^) can be Taylor 
expanded in Lp (or if we had applied the reduction of 
order procedure), to leading order, and for x ^ x , we 
would have obtained Q^, va p{x, x') ~ L p T : / lvaj3 (x — x')~ 4 , 
a result qualitatively similar to that of the analogous 
two-point function at one loop in perturbative quantum 



gravity. For solutions of semiclassical gravity with other 
scales present apart from the Planck scales, induced grav- 
itational fluctuations may occur on macroscopic correla- 
tion scales. In such cases, the above results suggest that 
stochastic semiclassical gravity might yield physically rel- 
evant results which cannot be obtained from a calculation 
in the framework of perturbative quantum gravity. 
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